Dynamics of Open Bosonic Quantum Systems in Coherent State Representation 
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We consider the problem of decoherence and relaxation of open bosonic quantum systems from a 
perspective alternative to the standard master equation or quantum trajectories approaches. Our 
method is based on the dynamics of expectation values of observables evaluated in a coherent 
state representation. We examine a model of a quantum nonlinear oscillator with a density-density 
interaction with a collection of environmental oscillators at finite temperature. We derive the exact 
solution for dynamics of observables and demonstrate a consistent perturbation approach. 
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I. INTRODUCTION 

When a small quantum system interacts with an en- 
vironment it undergoes irreversible processes such as re- 
laxation and decoherence. Environment-induced deco- 
herence is at the heart of the quantum-classical transi- 
tion [10. Classicality is an emergent property induced 
on subsystems by their environment. Under a variety 
of conditions, which are particularly easy to satisfy for 
macroscopic objects, decoherence leads to the selection 
of a small subset of quasi-classical states from within the 
huge Hilbert space. 

The dynamics of decoherence and relaxation in open 
quantum systems is usually studied with the help of the 
master equation for the reduced density matrix p of the 
system la , or by means of any of its possible unrav- 
cllings 3, including continuous quantum measurement 
and quantum trajectories g g & @ . It is also pos- 
sible to study the decoherence process in terms of phase 
space densities, especially used in the field of quantum 
optics |lfjj| . In this context one expresses the density 
operator p of a harmonic oscillator (one mode of the 
electromagnetic field) in terms of a c-number function 
of a coherent state (complex) variable a. Particularly 
useful distributions are: a) The Q-function, Q{a,a*) = 
— {a\p\a), that allows the calculation of antinormally or- 
dered quantum expectation values in terms of simple mo- 
ments of Q(a,a*); b) The P-function, P(a, a*), defined 
via p — J d 2 aP(a, a*)|a) (a|, whose moments give nor- 
mally ordered quantum expectation values; and c) The 
Wigner function W(a,a*), whose moments are equal to 
the expectation values of symmetrically ordered prod- 
ucts of creation and annihilation operators. These phase- 
space distributions have certain drawbacks. For example, 
they may no be positive definite, or may make no sense 
for certain density matrices. It may be difficult to ex- 
tract physical information from these quasi-probability 
distributions, especially in the context of nonlinear open 
quantum systems. 

The aim of this paper is to introduce an alternative 
approach based on c-number dynamical equations for ex- 
pectation values of observables of open bosonic quantum 



systems. The method is a generalization of the well stud- 
ied asymptotic theory for bosonic closed quantum sys- 
tems to the case of open dynamics. Such an approach 
provides a method to derive exact, c-number, partial dif- 
ferential equations describing the evolution of quantum 
averages evaluated in coherent states. In this sense the 
method is related to the phase space distributions dis- 
cussed in the above paragraph, but has the key feature 
of dealing directly with expectation values rather than 
with the quantum state p. This change of focus has sev- 
eral advantages. First, the physical interpretation of the 
results is clearer, and one does not need to compute inte- 
grals over moments to obtain physical quantities, which 
is especially difficult in the quasi-classical regime of pa- 
rameters due to fast oscillations. Second, the differential 
equations are well behaved in the quasi-classical limit 
h/N — > 0, where N is a quasi-classical parameter of the 
system, and lead to asymptotic Laplace-type expansions 
[l2| . The crucial property of Laplace asymptotics is that 
observable quantities arc exponentially localized in phase 
space around coherent states, and do not have the stan- 
dard oscillatory WKB behavior 01 . Third, it provides 
a simpler interpretation of singularities (in the sense of 
perturbation theory of partial differential equations [Tl| ) 
that appear in the quasi-classical regime for nonlinear 
Hamiltonians, allowing one to distinguish which part of 
the singularity is connected with measurable physical 
phenomena, and which is connected just with the choice 
of representation. 

As a prototypical system we consider the dynamics of 
a quantum nonlinear oscillator (QNO) 



Hs = fkjja'a + ph (by a) , 



(1) 



interacting with a bath of linear oscillators initially in 
thermal equilibrium. Here a (a^) arc annihilation (cre- 
ation) bosonic operators, to is the linear frequency and 
p is the parameter of nonlinearity. The QNO is initially 
prepared in a coherent state |a) in the quasi-classical 
region of parameters. The quasi-classical parameter is 
hj J <C 1, where J = fi|a| 2 is the action of the linear clas- 
sical oscillator. The nonlinear harmonic oscillator Hamil- 
tonian Eq. Q may describe a Bose-Einstein condensate 
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(BEC) treated in the single mode approximation. Such 
an approximation is valid when the many-body interac- 
tions within the condensate produce a small modification 
of the ground state of the trap, the mode structure be- 
ing sparce, such as in tightly optically trapped systems. 
For BECs trapped in optical lattices, the single-mode 
approximation leads to predictions for the quantum dy- 
namics of the condensate in excellent agreement with ex- 
periments 0,0]. Other systems that can be described 
by this Hamiltonian are micro- |T^ | and nanomechani- 
cal fl^ ] resonators in the nonlinear regime, and nonlinear 
optical systems, among others. 

The quantum and the classical dynamics of the an- 
harmonic oscillator was studied in detail in using 
the Q phase-space distribution of the system. It was 
shown that the presence of non-positive-definite second- 
order terms in the quantum evolution of Q, not present 
in the evolution of the classical probability distribution, 
is responsible for quantum recurrences and prevents the 
appearance of fine-scale-structure "whorls" predicted in 
the classical description. In [20I l2lj the interaction of 
the nonlinear oscillator with an environment (modeled 
by a thermal bath of harmonic oscillators with position- 
position coupling to the system oscillator) was studied 
in the limit of small nonlinearity using the Q function 
formalism, and it was argued that such an interaction 
was effective in destroying quantum interference effects 
and restoring the classical phase-space structure. How- 
ever, as recently shown by some of us in |l8l | by means 
of solving the master equation for the reduced density 
matrix of the system, environment-induced decoherence 
is in fact ineffective in recovering the quantum-classical 
correspondence for this nonlinear system: Some quantum 
effects may survive the decoherence process, and be ob- 
served for times much larger than the decoherence time- 
scale. In particular, we showed that the Ehrcnfcst time 
t E = (27tyt|a|) _1 , which characterizes the departure of 
quantum dynamics for observables from the correspond- 
ing classical dynamics, can be observed for times much 
longer than the decoherence time-scale. 

The paper is organized as follows: In Section II we con- 
sider a model of phase decoherence that is exactly solv- 
able, and will be used to exemplify our method based on 
the coherent state representation for observable values of 
open bosonic quantum systems, that we describe in Sec- 
tion III. This section contains the main results of this 
paper: The general theory proposed in this paper is first 
described, and it is then applied to the case of the non- 
linear oscillator, presenting both exact and perturbative 
treatments. Finally, Section IV contains our conclusions 
and briefly discusses possible extensions of this work. 



II. MODEL FOR PHASE DECOHERENCE IN 
THE NONLINEAR OSCILLATOR 

In this section we consider the model of phase decoher- 
ence in the nonlinear oscillator, that allows for an exact 



solution for the reduced dynamics of the system. This 
model will serve us for ease of presentation of the general 
theory of dynamics of open bosonic quantum systems in 
coherent state representation that will be described in 
the next Section. The interaction between the nonlinear 
oscillator and the thermal bath of harmonic oscillators 
with Hamiltonian 

N 

H E = J2^Mh, (2) 

3=1 

is of the density-density type 

jv 

H^t = ^2gja)ab]bj. (3) 

3 

Here bj (bij) are annihilation/creation operators of the 
environment, uij are the frequencies of the environmental 
harmonic oscillators, gj are coupling constants, and N 
is the total number of the oscillators in the environment. 
This model of decoherence can be used as an approximate 
description for the collisional effects between an ultracold 
atomic gas and its thermal cloud [2*3. f23| . Indeed, for low 
enough temperatures phase decoherence (corresponding 
to density-density type of interactions, i.e., elastic two- 
body collisions that conserve the total number of con- 
densed particles) dominate over amplitude decoherence 
processes (that correspond to position-position coupling, 
i.e., inelastic two-body processes that either feed or de- 
plete the condensate). 

The joint system-environment dynamics can be exactly 
solved in the Fock basis. Let us assume that the initial 
joint state is uncorrelated, i.e., /Otot(0) = p(0) <& Pb(0). 
The initial state of the system is assumed to be a pure 
coherent state p(0) = |a)(a|, and the initial state of 
the environment is assumed to be a thermal state p E = 
Z^ 1 e~^ B / kBT , where Z E = TT E [e^^ E/kBT ] is the parti- 
tion function of the environment. Here T is the tempera- 
ture of the bath, and fee is Botzmann constant. The joint 
density matrix at time t can be easily computed in the 
number representation since both the system, the envi- 
ronment, and the interaction Hamiltonians are diagonal 
in that basis. The result is 

oo 

Ptot(t)= Y, a n a* n ,e- u ^ n - n '^^ n2 - n ' 2 V\n){n'\® 

n,n'= 
N oo h (j) 

Z^H jr e -^e-^("-"V^O) )(/i G-)|. 

Here / Vnl are the expansion coefficients 

of the coherent state \a) in the Fock basis. 

Tracing over the environmental degrees of freedom it 
is easy to find the reduced density matrix for the non- 
linear oscillator p(t), and calculate any expectation value 
of a system operator. In particular, the evolution of the 
coherent state amplitude is given by 

(6(t)) - a(t)R(t), (4) 
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where a(t) is the solution without coupling to the envi- 
ronment, namely 

a(t) = a e -*(">+^)* e^[\a\ 2 (e- 2 ^ ht - 1)], (5) 

and R(t) = Ylj Rj(t) is the decoherence factor, that con- 
tains the effects due to the environment. Each Rj can 
be written in terms of its modulus and phase, Rj(t) = 
|iJ,-(i)|e*^W, where 



\Rj(t)\ 



tan ifj (t) 



1 - e fc B^ 



1 - 2e fc B^ cos(gjt/H) + e 
e fc B T sin(g,,-i/?i) 



1 — e fc B T cos(gjt/H) 



(G) 



The coherent amplitude a(t) departs from the classical 
solution on the Ehrenfest time scale ts = l/2Hfi\a\, and 
undergoes collapses and revivals, the revival time being 
The effect of the environment is to pro- 
duce a dephasing of the coherent state amplitude a(t), 
causing it to decay to zero in the limit of a large envi- 
ronment (N ^> 1), and killing the revivals. For finite N, 
the decay is incomplete and the revivals are suppressed. 
In the special case of identical environmental oscillators 
(gj = g and ujj = lue for all j), the decoherence factor 
R(t) can be approximated in the limit of large N by a 
periodic Gaussian structure. Each Gaussian R p (t) is cen- 
tered around a time t p = 2irhp/ g (p = 0, 1, 2, . . .) and has 
the form 



Rp(t) = cxp 



g 2 N{t-t p f 
2H 2 



(1-e fc BT)5 



III. COHERENT STATE REPRESENTATION 

There are several ways of considering the dynamical 
behavior of quantum expectation values. One possibility 
is to solve the Hciscnbcrg equation for the density matrix 
in some basis states and then take expectation values of 
the corresponding operators, as done in the previous Sec- 
tion. For general nonlinear Hamiltonians this approach 
leads, in the quasiclassical asymptotic limit, to a singu- 
lar behavior (in the sense of perturbation theory of par- 
tial differential equation) [Tj. Instead, it is possible to 
write down exact, c- number partial differential equations 
for expectation values, which allows one to distinguish 
which part of this singularity is connected with measur- 
able physical phenomena, and which is connected just 
with the choice of representation. This method is de- 
scribed extensively in pd| for closed bosonic quantum 
systems, and it is based on computing quantum observ- 
ables in the coherent state basis. In the following we 
briefly review the methodology for this case, and then 
we generalize it to open systems, i.e., quantum systems 
in interaction with an external environment. 



A. Closed systems 

In this subsection we exemplify the methodology of 
coherent state representation for closed systems using 
the nonlinear oscillator described with the Hamiltonian 
Eq.Q). Similar ideas can be applied to any quantum 
boson and spin systems, as described in Given an 

arbitrary operator of the system fs = fs(a\a), it is 
possible to write down an exact, c-number partial differ- 
ential equation for the time-dependent expectation value 
fs(a*,a,t) = (a\e iAst / R f s e- i6st / h \a) of such operator 
evaluated in coherent states |ct). Using Heisenberg equa- 
tion d/dtfs — (i/h)[Hs, fs], it follows 



fs = \ ((a\H s f s \a) - (a\f s H s \a)). 



(7) 



Putting both Hs and fs into normal-ordering form in 
terms of the initial operators a = a(t = 0) and = 
d^(< = 0), it is possible to obtain a closed form for each 
of the two terms in Eq.(0. One obtains 



(a\H s fs\a) 
(a\f s H s \a) 



H s a 



4 2 



fs e 



M 2 



(8) 



Here fs = fs(t) is the c-numbcr expectation value we 
seek, and the differential operator Hs has the same 
functional form of the normal-ordered operator func- 
tion H 5(0,^,0,), but with the substitution a) — > a*, and 
a — > Therefore, the exact partial differential equa- 

tion for the time-dependent expectation value fs(t) reads 



- jrr = K s fs 



dfs 
dt 

with initial value /s(0) = fs(& 
ferential operator Ks given by 

d 



(9) 

and the partial dif- 



K f 



H s a 



ck 



-H s [a 



d_ 

da 



This operator can be split into two parts: Ks 
HK q . The first operator K c \ includes only 



>l . 

(10) 

= Kci + 
first or- 
der derivatives and describes the corresponding classical 
limit, while the second operator K q includes higher or- 
der derivatives and is responsible for quantum effects. 
For the model described by Eq. (JTJ the exact partial 
differential equation for observables reads 



dfs 
dt 



= i(oj + h[i + 2ft/i|a| 2 ) I a 



d 2 



ihfi (a*Y 



d(a*f 



,_d 
da 

da 2 







fs + 



fs. 



(11) 



In particular, for fs = a, the evolution of fs(t) cor- 
responds to the evolution of the coherent amplitude 
a(t) = (a\a(t)\a), and the solution is the same as in 
Eq. ©. 
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B. Open systems: Exact treatment 

In order to extend the formalism to treat open quan- 
tum systems, we assume that the system, initially popu- 
lated in a coherent state \a), interacts with the bath of 
harmonic oscillators, also initially populated in coherent 
states {l/Sj)}- Let / be any operator of the composite 
system 



/ = /(at (i ),a(t), {£].(*)}, &.(*)}), 



(12) 



that evolves according to the Hcisenbcrg equation 
df/dt = | [H,f], where H = H s + H E + H int is the 
Hamiltonian of the composite system. The initial state 
is |^(0)) = \a,{f3j}). The expectation value f(t) of the 
operator f(t) in the state |^(0)), evolves according to the 
partial differential equation 

df 
dt 



(13) 



where K is the differential operator 



K = -e^-^^M 

h 



d_ 

da ' 



2 +E 3 Iftl 2 



(14) 



Note that the following substitutions have been used: 
at -> a*, a gfy, b) [3*, and bj -> For 

the model under consideration, the c-number differential 
equation for the composite (system+bath) expectation 
value f(t) has the form 



df 

= {K a +K, 3 + K int )f. 



(15) 



The first two terms in Eq. I|15|) correspond to the free 
system-bath dynamics, 



K a f = i{u + nh + 2pLh\a\ l ) [ a*^--a-^- ) f 



ifj.h 



d 2 



d(a*) 



*\2 



da 
da 2 



da 



(16) 
(17) 



and the last term in Eq. (|15fl is due to the system-bath 
interaction. It can be written as the sum of three contri- 
butions K int f = (Kg* + Kg> + K$>)f, where 



k {1) f 

ly nt J 



-"-int J 



-"-int J 




|2 \ " 



d 



a 



0_ 

da 




3 



fJj 



* d d 
da* Pj d(3* 



a-^(3 



da ' 3 d/3j 



Given a solution /(a, a*, {f3j}, {P*}\ t) to Eq. {TBJ), we 
finally have to trace over the coherent states 
(i.e., trace over the environment) to obtain the evolution 
of expectation values of the system. 

In the above we have assumed that initially each j-th 
environmental oscillator is in a pure coherent state \(3j). 
Let us now consider the case in which each environmental 
oscillator is initially in a mixed, thermal state at temper- 
ature T. Since the oscillators in the environment are non- 
interacting, the initial density matrix of the environment 
can be written as a direct product over individual density 
matrices for each j-th sub- environment, that is p_e(0) = 

IljLi Pe ' (0)- The density matrix of each environmental 

oscillator is then p^(0) = (Z^)- 1 exp ( — hwj b^bj /ksT), 

where zf = Tr[e"^ & ^ j/fcl3T ] is the partition function 
of the j-th environmental oscillator. This mixed thermal 
state pjp(0) can be written in the coherent state basis 
(this corresponds to the so-called P-representation or co- 
herent state representation [24[) 



^(0)- / d^jPiP*,^)]^)^], 



(18) 



where the probability distribution for each & is given by 



P(j3*,/3 j ) = ^e 
J 7m,,- 



1 -'ftl 2 /^ 



(19) 



Herefij = ( e Hu; ^ kBT -l)' 1 i s the Bose distribution. That 
is, the coherent state representation of a thermal state 
has a Gaussian distribution. 

For the particular case of the system oscillator initially 
prepared in a coherent state, we would like the solution 
to Eq. (|15fl to lead us to Eq. once the integration 
over the environmental variables is performed. Since the 
structure of Eq. is the product of the free evolution 
solution a(t) times a time-dependent factor that arises 
from the coupling with the environment, we propose a 
solution to Ea. l|15|) of the form 

/(a, a*, {/?*}; t) = f a (a, a*;t)x />({/%}, {/?*}; t), 

(20) 

where f a = a(t) is given by Eq. JSJ). The reduced dy- 
namics of an observable of the system will be given by 
an average over the environmental oscillators weighted 
by their respective probability distributions 



fs{t)=a(t) x 



The initial condition for fp(t) is fp(t = 0) = 1. Since 
P(/3*,/3j) depends on (3j through its modulus squared, 
we can assume, without loss of generality, that at any 
time t the function fp(t) depends on \/3j\ 2 (any other 
dependency, like f3j((3*) m (k ^ m), vanishes upon inte- 
gration). This implies that the operator Kp acts on fp 
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as Kpfp = 0. Given that f a = K a fa, the equation for 
fp finally reads 



fp = fa Kixttfafp- 



(22) 



As can be shown by direct inspection, an exact solution 
to this equation, with initial condition fp(t = 0) = 1, is 



//j(*)=ntf ) (*)=ii e ^H^i a ( i - e 



igjt/h 



(23) 

Integrating upon the probability distributions P{\Pj\ 2 ) 
we obtain the reduced dynamics for the observable of 
the system oscillator, namely 



f s (t)=a(t)l 



d 2 P P{\P 3 \ 2 )f? 



a(t)R(t), (24) 



where R(t) = . Rj if) is given in Eq. © . This expres- 
sion coincides with the exact solution derived in Eq. 
obtained from solving the Heisenberg equation for the 
density matrix of the joint quantum nonlinear oscillator- 
environment system. Therefore, using our method based 
on the coherent state representation for observables val- 
ues for open bosonic quantum systems, we can recover 
the exact reduced dynamics. 



do not have exact solutions, we will now solve Eq. i|15|) 
for this model considered in a perturbative expansion in 
powers of H lrA , and compare the results with the exact 
solution previously found. The study of other nonlinear 
models will be left for a future publication. 

The perturbative master equation for the reduced den- 
sity matrix of the quantum nonlinear oscillator is, to sec- 
ond order in H mt , 



dt 



p = — i[(ui + 8u>)n + phh 2 , p] + tj(2npn — h 2 p — pn 2 ), 



where h = a^a is the number operator for the sys- 
tem. The first term is the free unitary evolution 
with an environment-renormalizcd frequency 5lo = 
(1/h) J2j 9j{nj}> that arises from first order perturbation 
theory. The second term is of Lindblad form (but with 
a time-dependent coefficient), it arises from second order 
perturbation theory, and it is responsible for decoherence. 
The coefficient 7 is given by 7 = (1/h 2 ) J^. g 2 [(n 2 ) — 

(fij) 2 ], where fij is the number operator for the jth oscil- 
lator in the environment. The perturbative treatment in 
powers of H mt is in fact a perturbative treatment valid 
for short times {t/h)(n)Y^,j Qjftj) ^ !• The solution 
of this approximate master equation can be straightfor- 
wardly found in the Fock basis, and from there one can 
evaluate the dynamics for the coherent state amplitude 



C. Open systems: Perturbative treatment 

There are a few models of decoherence and relaxation 
for which it is possible to solve exactly the joint system- 
environment dynamics, and to write down and solve an 
exact master equation for the reduced density matrix of 
the system. Examples are the one considered in Section 
II, that can be trivially solved in the number represen- 
tation, and the well-known quantum Brownian motion 
model, in which a linear oscillator (p = in Eq.QJ) is 
coupled through position to a bath of linear harmonic 
oscillators. This latter model can be solved using, for 
example, influence functional techniques j2^| thanks to 
the fact that both the system, the bath, and the interac- 
tion Hamiltonians are quadratic forms. For other general 
models of decoherence and relaxation, and in particular 
for nonlinear open quantum systems (for example, the 
nonlinear oscillator of Eq. (JJJ coupled via position to the 
bath of linear oscillators), there are no known exact solu- 
tions. In those cases it is customary to use different ap- 
proximation methods, such as a perturbative expansion 
in powers of the interaction Hamiltonian H mt (Born ap- 
proximation), and, when applicable, the Markovian ap- 
proximation (memoryless environment) [Tol |. 

The model considered in the previous sections (the 
nonlinear oscillator coupled via density with the bath 
of linear oscillators) affords an exact solution due to the 
simplifying property that all terms in the Hamiltonian 
are diagonal in the joint (system+bath) number basis. 
In order to show how to deal with exact PDEs for ob- 
servables in generic bosonic open quantum systems that 



(a(t)) 



a{t). 



(25) 



We see that the coherent state amplitude has an initial 
quadratic time decay, typical of quantum systems sub- 
jected to perturbations. As expected, this perturbative 
solution is the same as that obtained from the exact so- 
lution presented in Eq. Q. Indeed, for short times the 
modulus and amplitude of each decoherence factor Rj (t) 
can be written as 



\Rj{t)\ 



1 



g*t 2 



2h 2 



[in 2 ,) 



9it 



(26) 



Therefore, the total decoherence coefficient is given by 



R(t) « [1 - 7* 2 /2]e 



-iSut 



, that coincides 



with the solution of the perturbative master equation. 

We now turn to find the perturbative solution to the 
exact PDE for observable values, Eq. l|15|) . Inspired in 
the solution above, we again propose a solution of the 
form f(t) = f a (t) fp{t), with f a (t) = a(t). Let us take 
a perturbative expansion of fp(t) in powers of the small 
parameter e of the form 



AO) 



Jp 



(27) 



where is independent of e, /g 



Given the initial condition fp(t 
The first order equation is 



is linear in e, etc. 

•(0) _ 1 

P ~ 1 - 



0) = 1, then A 



Jp 



fa K'mtfctf, 



(0) 



(28) 
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whose solution reads / 
to first order we obtain 



(i) 





-it 



Ejftl&f- Therefore, 



it 



(29) 



Integrating upon the probability distributions P((3*,f3j) 
we recover the short-time solution Eq. (|25[) . 

fs(t) = a(t)x J]l*0iP(fi,Pj)fft{t) 



e 2 e 



-iSujt 



a{t). 



(30) 



which, again, coincides with the perturbative solution 
Eq. 

For other models of decohcrence, such as high temper- 
ature quantum Brownian motion, one can proceed along 
similar lines, i.e. solve the exact PDE for observables 
in a perturbative expansion in powers of the interaction 
Hamiltonian -Hint- It is possible to introduce an infinitely 
divisible distribution Y whose probability distribution is 
not Gaussian, but given by a Levy distribution, that leads 
to different time dependencies of the decoherence factor. 
For example, for a Lorentzian probability distribution 
one obtains an exponential decay [27J |28|, |2!j . 



Another way of obtaining the same result is to use 
concepts of probability theory. This may turn out to 
be useful in other models of decoherence for which long- 
time solutions for the reduced dynamics of the system 
are available [2(|. Let us call Xj = gj\Pj\ 2 the stochastic 
variable that takes the values Xj = gjUj with probability 

where fij is given by the Bose 
g 3 n = gj(xi), 
fti 



P{n ) = (nn^e-^/* 



distribution. The mean value of Xj is Oj 



and its variance is b 2 = (x 2 ) — (x,) 2 = g 2 [{^ 



— Vj/ — ~ ujW'ji ~~ \ n 'ji j- 
The perturbative solution fp(t) can then be written in 
terms of these stochastic variables as 



(31) 



The stochastic variables Xj can be considered as indepen- 
dent and identically distributed. Therefore, the stochas- 
tic variable Y = YljLi x j belongs to the class of the so- 
called infinitely divisible distributions [H HI H| . The 
behavior of Y depends on whether the cumulative vari- 



ance Bjj 



J2j bj is finite or not. In the limit N 



B N 



CO, 



is finite (central limit theorem), and the probability 
for Y is Gaussian 



P(y) = 



i 



exp 



{y - vf 

2B N 



(32) 



where y = J2j x j = J2j 9j( n j) ■ To obtain the reduced 
dynamics for the system we need to integrate fp(t) over 
the environmental variables f3j weighted with their prob- 
ability distributions. This is equivalent to integrating 
e -ity/h over j£ g probability distribution P(y), 



IV. CONCLUSIONS 



In this paper we have generalized the method of ex- 
act partial differential equations for observable values of 
bosonic systems 0, i based on a coherent state rep- 
resentation, to the case when the system interacts with 
a bosonic environment. Our method requires to solve, 
either exactly or approximately, a PDE containing cou- 
pled coherent state degrees of freedom of the system and 
the environment, and then to integrate (trace) over the 
environmental coherent states weighted by their respec- 
tive probability distributions. We have exemplified the 
method with a model of a nonlinear oscillator interacting 
via density with a bath of linear oscillators. The simplic- 
ity of this model, based on the fact that all term in the 
Hamiltonian arc simultaneously diagonal in the number 
basis, allows for an exact solution. We demonstrated that 
the dynamical behavior obtained from the coherent state 
representation coincides with that obtained from the re- 
duced density matrix approach. Further development of 
our method is required in order to treat other more com- 
plicated decoherence models, such as a nonlinear oscilla- 
tor coupled through position to the environment. This 
will involve the study of consistent perturbative solutions 
to the exact PDE for observables. 



dy e-**» P(y) 



(33) 
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fs(t) 



e- lSuJt a(t), 



(34) 
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